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We revisit the Hamiltonian formalism for a massive scalar field and study the particle 
production in a de Sitter space. In the invariant-operator picture the time-dependent 
annihilation and creation operators are constructed in terms of a complex solution to 
the classical equation of motion for the field and the Gaussian wave function for each 
Fourier mode is found which is an exact solution to the Schrodinger equation. The 
(N- in-out formalism is reformulated by the annihilation and creation operators and the 

Gaussian wave functions. The de Sitter radiation from the in-out formalism differs from 
the Gibbons- Hawking radiation in the planar coordinates, and we discuss the discrepancy 
■ of the particle production by the two methods. 
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■ 1. Introduction 



The universe with a cosmological constant is the pure de Sitter (dS) space and 
has the maximal symmetry that makes a quantum field theory more tractable than 
other curved spacetimes (for a review and references, see ref. [I]). In spite of numer- 
ous works since the discovery of the dS-invariant vacuum by Bunch and DavieJ^l 
and the dS radiation by Gibbons and Hawking (for instance, see ref. [4]), the 
proper selection of the physical vacuum of a quantum field in dS space has recently 
been challenged by Polyakov.^^ It has been further argued that the cosmological 
constant problem may be resolved by IR physics through the UV/IR mixing from 
an mteraction.'ii' 

In this paper, we revisit the quantum field theory for a massive sc alar field in 
the dS space by unifying the in-out formalism by Schwinger and DeWitlP^^ and the 
invariant-operator picture by Lewis and Riesenfeld.^^ The in-out formalism provides 
a good framework for quantum field theory involving particle production. It has 
been known that an expanding spacetir ne p roduces particles-'^'^ and that the dS 
space has the Gibbons-Hawking radiatiorP'^ (for a recent discussion, see ref. [11]). 



1 



April 19, 2012 0:23 WSPC/INSTRUCTION FILE cospa2011-spkini 



2 S. P. Kim 



Recently t he dS rad iation has been interpreted in connection with the Schwinger 
mechanism! ^1^1^ 2 113j contrast to the early works on the in-out formalism jl3! in this 
paper we shall directly use the wave functions for the field and shall find the vacuum 
persistence amplitude between two spacelike hypersurfaces, one in the past infinity 
and the other at an arbitrary time. In cosmological scenarios the Gaussian wave 
functions on each hypersurface carry the same information as the Heisenberg field 
operator. In fact, the in-out formalism and the invariant-operator picture formalism 
are equivalent to each other. 

In the Hamiltonian formalism a massive scalar field in the dS space is equiv- 
alent to an infinite sum of time-dependent harmonic oscillators.^^ The invariant- 
operator picture provides a pair of time-dependent annihilation and creation op- 
erators t hat generate all t he exact quantum states in analogy with a harmonic 
oscillator! ^^ l ^^ l ^^ l ^^ | 20 | 2 1] Uging these operators we construct the Gaussian wave 
packets for each oscillator that evo lve from the sp acelike hypersurface in the past 
infinity toward to the future infinity. P'^ l '*^^ l -*^^ l 20 | 21| p^-|-|-]-^g-,.j^Q-|.g^ .^^g gj^^j ^t^q Bogoli- 

ubov transformation between the operators on the two different hypersurfaces. 

Interestingly, in the planar coordinates of dS space the number of in-vacuum 
particles carried by the Gaussian wave packet on each hypersurface exponentially 
increases in proportional to the cube of the scale factor and the vacuum persistence 
amplitude confirms this phenomenon. This contrasts with the Gibbons-Hawking 
radiation which is a thermal and scale-factor free distribution in the future infin- 
ity. It is shown that the matrix amplitude between the Gaussian wave packet and 
the out-vacuum Gaussian wave packet in the future infinity, however, exhibits the 
Gibbons-Hawking radiation. The Gibbons-Hawking radiation is the particle number 
measured by any detector whose quantum states are prescribed by the out-vacuum 
solution. Hence, from the view of the in-out formalism this discrepancy raises a 
fundamental question on the quantum field in the planar coordinates of the dS 
space. 

2. Hamiltonian Formalism for Massive Scalar Field 

In a Friedmann-Robertson- Walker spacetime [in the units of c = ?i = 1] 



ds^ = dt^ 



a'^{t)dx^, 



(1) 



a massive scalar field $ has the action 





In the Hamiltonian formalism the action with H = 0^^$ 




(3) 



is given by the Hamiltonian and its density 
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The stress-energy-momentum tensor 

T^, = = a^<i>9.$ - ig^p (^gP'-dp'S^d,^ ~ m^*^) , (5) 

has also the Hamiltonian expression 

Too = ^, To^ = 

Now we decompose the Hamiltonian (U) by the Fourier expansion of the field 
and the conjugate momentum 

n(t,x) . j j^,Mt)e^-- (7) 
Then the Hamiltonian is an infinite sum of time-dependent oscillators 

Hit) = J d'k^[^ + a^u^Ul) = J d'kHkit), (8) 

where 

TTk = a^(t)k, ujl{t) = — +nt'. (9) 
Here and hereafter (/>| and 7r| denote 

: = -k.r_k = j + j ' (10) 

and fc2 = k^. Thus, the fc-mode denotes both the cosine and the sine modes. The 
Hamiltonian (|5]) has the meaning of Hk{t) = HkiT^k, 4>k, Sf) on each spacelike hyper- 
surface Sf and the Hamiltonian formalism describes the evolution from one hyper- 
surface to another Sf . Similarly, the stress-energy- momentum tensor ([B]) can be 
decomposed into the Fourier modes according to eq. (1101) . The Hamilton equation 
for Hk (t) is the Fourier mode of the field equation 

4>k + 3-(j)k+cvl{t)(t)k =0- (11) 
a 

3. de Sitter Radiation of Massive Particles in Planar Coordinates 

The maximal symmetry of a dS space allows one to explicitly solve the field equation 
for a massive or massless scalar field. In fact, in the planar coordinates, a{t) — e^*, 



ds^ = ~df + e^"*d:>i^, 



(12) 
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eq. pT|) in the canonical form takes the form of the Morse potential while in the 
global coordinates, a{t) — cosh{Ht) / H , its is related to the Posch- Teller potential. 
The general solution to eq. pT|) in the metric (IT2t is given by 



\4:HJ 



We--/^<(z) + cfe-<(z) 



(13) 



where H^p' and H^^'' are the Hankel functioni^ and for massive particles (m/H > 
3/2) 



k fm^ 9\i/2 
z=—e"\ p=l 



The standard quantization rule requires that the Wronskian condition should be 
satisfiecPl 



' (i)Wr[^fe (t) , it)] = a\t) (^fe mi it) - {tWl (t)) 



(15) 



Though the integration constants c\^^ and c^^^ can be arbitrary complex numbers. 



(16) 



the quantization rule (|15p restricts them to satisfy 



(1)|2 



„(2)|2 



1. 



In the past infinity (t ~ — oo), the general solution (jl3l) has the asymptotic form 

^.(t)=4\f)(t)+ci^vr^*w, (17) 

where the in-vacuum (the Bunch-Davies vacuuni^ with c^, — 1 and cj, — 0) is 
provided by the positive frequency solution with respect to the energy operator, idt, 



exp I — e 
2k H 



(4-"*)- 



(18) 



In the f utu re infinity (t = oo), using the asymptotic formulae for the Hankel 
function^ 



1 r e-^p 



sinh(7r]3) Lr(l + ip)\2 



1 



r(i -ip) \2 



H\l\z) = J,p{z)-zN,p{z) 



1 



z\v 



1 



z\ ~w 



■(I) 



(19) 



sinh(7rp) Lr(l + ip) \2J ' r(l - ip) 
the solution (fT3|) asymptotically becomes 

Mt) = awt'Ht)+p,^^r'^*{t). (20) 

where the out-vacuum is constructed by another positive frequency solution 

g-3Ht/2 



(out) / ,s 



'iHpt 



V^Hp 



(21) 
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Here the Bogoliubov coefficients are given by 

sinh(7rj3)r(l + ip) ' 

//VJ1)„-^p/2 , ..Jl)..p/2^ (fc/2g)-'P 

sinh(7rp)r(l — ip) 
A direct calculation shows that the Bogoliubov relation holds 

Kl^ - 1/3.1^ =(|4^^|^- 14:^^1^)=!. (23) 
Noting that is the particle number, the probability for particle production 

is 

^^■ = ^2^l4^^P + ^2i^|cfr (24) 

In the above we have assumed the real c[,^' and c^,'^' for simplicity which corresponds 
to the zero squeezing angle. The interpretation of eq. (I24p is that the first term is 
the spontaneous production while the second term is the stimulated emission in the 
presence of in-vacuum particles. In the specific case of the Bunch-Davies vacuurn, 
the probability for particle production is the Gibbons-Hawking radiation of bosona^l 

- ^5^- (25) 

The particle production is independent of the momentum and the scale factor a(t) 
and is determined by the mass and the Hubble constant. 

4. Gaussian Wave Packets 

We shall use the Schrodinger picture since quantum states are c-number wave func- 
tions that carry the same information as the Heisenberg field operator. Though not 
considered in this paper, the Schrodinger picture is also convenient for a nonlinear 
theory of self-interactions, modulo the renormalization of the wave function, the 
vacuum energy and the coupling constants. 

The time-dependent Hamiltonian (|8]) in an expanding universe may give an 
adverse feeling for the Schrodinger picture 

j^«'fe(</'fe,i) = Hfe(t)«'fe(0fc,i). (26) 

However, it has been known for a long time that the Gaussian wave function could 
be expressed in terms of a classical solution to eq. (fTT|) . for instance, see ref. [23] . 
It may be understood from the fact that the Heisenberg field operator also satisfies 
the same equation (jlip and the Heisenberg picture provides the same quantum 
information as the Schrodinger picture. In fact, the Gaussian wave function on a 
spacelike hypersurface of constant t for each Fourier mode, which is the solution 
to the Schrodinger (|26p and is normalized to unity, is given bjff^ 

(27) 
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where ipk is a complex solution to eq. (jlip that satisfies the quantization rule (jl5l) . 

In the invariant-operator picture by Lewis and Riesenfeld, the Hilbert space for 
a time-dependent Hamiltonian of quadratic order can be constructed exactly in the 
same manner as for a time-independent Hamiltonian. In this picture we look for the 
operators that satisfy the Liouville-von Neumann equation 

i^&feW + Mt)Hk{t)] = 0. (28) 

Note that eq. ((28)) is the equation for the density operator and describes the back- 
ward evolution 

b'^k\t) = ulit)bf^mt), (29) 
in contrary to the forward evolution of the Heisenberg operator 

U^\t)^U,{t)bf^ul{t), (30) 
where Uk is the evolution operator 

i^Ait) - Hk{t)lJk{t). (31) 
Indeed the operator satisfies eq. with respect to the Hamiltonian 

H^\t)^Ul{t)Hf^Mt)- (32) 

For the quadratic Hamiltonian we can directly find the invariant operators (the 
superscript (I) dropped be low') for 60- which are first order in the momentum 
and the position operatorsjlSllZIIHIIIIinill] 

bu{t) = i{ipl{t)%k - {t)ipl{t)li>k) , 

bl{t) = -i{ipk{t)frk - a^t)^k{t)^k). (33) 

In the above itk and 4>k are the Schrodinger operators and ipk (t) is a complex solution 
to eq. pT|) together with the quantization rule p^ . Remarkably, these operators 
satisfy the equal-time commutation relation 

[6fe(t),St(t)] = l, (34) 

and, moreover, play the role of the time-dependent annihilation and creation oper- 
ators since they reduce to the standard ones for a time-independent Hamiltonian, 
modulo time-dependent phase factors. In fact, the Gaussian wave function (P7| is 
annihilated by bk{t). 

5. Quantum Dynamics 

The Hamiltonian ([8]) has the algebra SU(1,1), for which we may choose a time- 
independent Hermitian basis 

L^-'^l4>l L(+^^hl L(o) = l(0,fr,+^fe<^,). (35) 
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As the field and the conjugate momentum operators have the expression 

k{t) = ^k{t)hk{t) + ip*k{t)bk{t), %k{t) = a^{t) [ipk{t)bk{t) + ipUt)bk(t)] , (36) 

OA 

we may use another Hermitian basis 

Mi"\t)^bl{t)h{t) + bk{t)bl{t), (37) 

Mt\t) = blit)+bi'it), Mir\t) =^[bl{t) bWt)]. (38) 

This basis has the SU(1,1) group structure 

[M^^\t),Ml^\t)] = ±2iM^^\t)5kk', [Mi-^Ht),Mi7\t)] = -2zM["\t)Skk' ■ (39) 

Note that Mj^^\t) is the number operator on the hypersurface S*. 

Using eq. (155)) we express the quadratic operators for the field and the conjugate 
momentum in terms of the basis p7l [55]) 



1 



>(+)/ 



+ \{^f{t)-^l{t))Mi-\t) 



Conversely, the number operator can be expressed in the basis (1351) as 



(40) 



The quadratic variances with respect to the Gaussian wave function (|27p on take 
the values 



{^kml\^k{t)) - Wt{t)\\ {^^{t)\7:l\^,{t)) = a\t)\^,{t)\\ 
{^k{t)\^kT^k + vrfc0fe|*fe(t)) - a3(t)(|^fc(t)n-. 



(42) 



6. Particle Production in Quantum Dynamics 

We may raise the question whether the operator p7p indeed provides a proper 
definition for the particle number 



The Gaussian wave function (j27l) has the zero particle number on Et 

{^k{t)\Nk{t)\^k{t)) =0. 



(43) 
(44) 
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Thus the number operator 7V^ {t) measures the amount of particles carried by the 
Gaussian wave function on each Et through the evolution, which is a conserved 
quantum number, that is, the harmonic excitation number. We use the Gaussian 
wave function on the hypersurface S_oo hi the past infinity as the in-vacuum state 
for each k and define the in-vacuum as the tensor product 

|vac,in) =[]|*fe(-oo)). (45) 

k 

Similarly we define the out- vacuum as 

|vac,out) =]J|*fc(oo)). (46) 
fe 

Then the number of in-vacuum particles contained in the Gaussian wave function 
on Et or the number of out-vacuum particles contained in the in-vacuum Gaussian 
wave function on E_oo is 

Ukit) ^ (*fe(i)|7Vfc(-oo)|*fe(i)) = (^'fc(-oo)|iVfc(i)|«'fe(-oo)). (47) 

The second equality is a consequence of the reciprocality, as will be explicitly shown 
below. 

Taking the expectation value of the operator Mj,^\—oo) in eq. ((4T|) with respect 
to the Gaussian wave function on S(, we obtain 

nk{t) = a%t)\M-<x)\^\v^kit)\^ + a^(-oo)|<^fc(-c3o)n(^fe(t)|2 

-^a'i-oo)a'mM^ooryi\Vkit)n \. (48) 

We find the Bogoliubov transformation between and 

ak{to) = Hk{to,t)ak{t) + i^k{to,t)al{t), 

aUto)=^^Uto,t)alit) + ,^*kito,t)akit), (49) 

where 

fik{to,t) = i[a^{t)ipl{to)(pk{t) - a^{to)'Pk{tWk{to)], 

Mto,t) = t[a^t)^UtoMit) - aHtoWimiito)]. (50) 

The Bogoliubov transformation (HH)) is symmetric under the interchange of ak{to), 
a|(io) and ak{t), al.{t), so the inverse transformation 

ak{t) = fJ.k{t, to)afe(^o) + i^kit, to)al{to), 

4(0 = ^Jl{t,tQ)al{tn) + vl{t,tQ)ak{tQ), (51) 

has the coefficients 

i^k{t,h) = /ifc(io,0) Vk{t,tn) = ~Vk{to,i). (52) 
It can be shown that the particle number is given by 

nk{t) = |iyfc(-(X3,t)|2 = |iyfe(t,-c3o)|^ (53) 
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Furthermore, the Bogohubov relation holds 

Mo,t)\^ ~Wkito,t)\^ ^1. (54) 

However, the Bogohubov coefficients ([^^ cannot be obtained from those ([5U|) using 
the asymptotic solutions in section [3] 

7. Vacuum Persistence Amplitude 

One important concept for the particle production and the vacuum in stabi lity is 
the vacuum persistence amplitude introduced by Schwinger and DeWitt,!^!^ which 
is the scattering-matrix (S'-matrix) amplitude between two different hypersurfaces 

e''^''^'^'^^ = i^kimkito)). (55) 

In particle physics the S'-matrix takes two asymptotic limits to = —oo and t = oo, 
that is, before and after interactions. But it is still legitimate to use the defini- 
tion (j55p in an intermediate region, though the physical meaning requires further 
clarification and investigation. The integration of two Gaussian wave functions (|27p 
gives 

Wk{t, to) = i \n{^lk{to, t)) = i ln{^ll{t, to)) , (56) 

where fj,k is given in eqs. ((5(11 [5^ and the mode k counts both k and — k. 
The vacuum persistence amplitude for the field is the product of eq. (|55p 

e'"^^'''"^ ^Yli^kimkito)), (57) 

k 

and thereby the effective action is 

Wit,to) = iVit) J d''klll{^lkito,t)) = iVit) J d3kln(/i*(i,to)). (58) 

Here iJ-k{to,t) is the Bogohubov coefficient obtained by substituting k for k in eq. 
([50)1 and V{t) is a spatial volume relevant to the hypersurface on Sj and may be 
chosen as V(<) = a,'^{t) J d'^x. The effective action ([55]) involves the divergent terms 
that renormalize the vacuum energy or the mass, and the gravitational constant 
or the charge. In the in-out formalism equation eq. ([55]) has been used to find the 
effective action for the massive scalar field in the global geometry of dS spacJI^I 
and in quan tum electrodynamics (QED) the effective action in a time-dependent 
electric field.l^SI 

The vacuum persistence (twice of the imaginary part) relates to the particle 
number through the relation 

2lm{Wkit,to)) = ln{\^ikit,to)\^) =ln{l + \iykit,to)\^). (59) 
And the probability for the in-state on to remain in the out-state on T,t is 

|(0,t|0,to)P = e-2MM^(*'*°»=exp[-V(i) / d'kln{l + nkit,to)\^)]. (60) 
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Thus the vacuum of dS space becomes unstable due to the dS radiation. It is shown 
that the number of in-vacuum particles carried by the Gaussian wave function on 
at a later time hypersurface St is dominated by the first term in eq. (I48p 

Ukit) a'{t){a'{mk{t)\')\Mt)\'\M-^)\' (61) 

and is enhanced by the exponential scale factor and by the infinitely blueshifted 
wavelength on the initial surface while the parenthesis is finite. This catastrophic 
production of in-vacuum particles in the infinity future is a generic phenomenon in 
an expanding universe with a cosmological singularity. An interesting observation 
is that if a regularization scheme removes a'^{t) and \(pk{—oo)\'^ in the vacuum 
persistence probability (I60p . the remaining term is related to the Gibbons-Hawking 
radiation. 



8. Wave Packet for Out- Vacuum 

The out-vacuum solution, 

(out) 



1/2 

2H sinh(7rp) 

satisfies the quantization rule (jl5p and has the same time-dependent factor as eq. 
([20)1 up to a constant phase factor. The corresponding Gaussian wave function 

(out) T /o ■ • (out)* 

is another exact solution to the Schrodinger equation (j26l) since eq. (j62[) is a complex 
solution to eq. ([Tl]) . The Bunch-Davies vacuum solution ^ (c^^' = 1, = 0) 
now has the Bogoliubov transformation 

Vk{t) = afc^[™*Hi) + /3fc^i°"*^* (64) 

where 

"'==i2sinh(.p)J ' ^'= = -l2^h^hMJ • 
The particle production |/3/tp is the Gibbons-Hawking radiation for massive bosons 

(ESI. 

In analogy with the vacuum persistence amplitude between the in-vacuum and 
the out-vacuum in section [71 we may compute the matrix amplitude on the same 
hypersurface St 

e'<""(*) ^(M,,(t)|vl;(™')(i)). (66) 

Then the vacuum persistence (twice of the imaginary part) relates to the particle 
number through the relation 

2Im(V7i°"*)(t)) = ln(|afe|2) = ln(l + (67) 

and explains the Gibbons-Hawking radiation. 
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9. Conclusion 

In this paper we have studied the production of massive particles in a dS space 
using the quantum dynamics that unifies the invariant-operator picture with the 
in-out formalism. In the quantum dynamics a massive field in an expanding FRW 
universe or a dS space is equivalent to an infinite sum of time-dependent oscillators 
in which the time-dependent mass comes from the expanding spatial volume and 
the frequencies from the redshifted or blueshifted wavelengths. The main advantage 
of the invariant-operator picture is that there exist the time-dependent annihila- 
tion and creation operators which generate the exact quantum states, the Gaussian 
wave function being the simplest one among them. The Gaussian wave function 
provides the in-vacuum in the past infinity and the out-vacuum in the future in- 
finity. Hence the time-dependent annihilation and creation operators connect these 
two asymptotic regions through the Bogoliubov transformation. 

We have found the Gaussian wave packets corresponding to the Bunch-Davies 
vacuum and its one-parameter family. The scattering-matrix of Gaussian wave pack- 
ets between two different spacelike hypersurfaces directly gives the vacuum persis- 
tence amplitude, which is equivalent to that from the Bogoliubov transformation 
and its coefficients. The vacuum persistence amplitude is an important tool in quan- 
tifying how the massive field probes the background dS space in analogy with QED 
in a strong electric field background. However, in the planar coordinates of dS space 
the number of in-vacuum particles carried by the Gaussian wave packet exponen- 
tially increases in proportion to the expanding spatial volume, in contrast with the 
Gibbons-Hawking radiation. This is because the momentum variance, the leading 
term for particle number, exponentially increases in the infinity future while the 
wave packet is very sharply peaked in the field. On the other hand, the matrix 
amplitude between the Gaussian wave packet and the out-vacuum Gaussian wave 
packet on the same hypersurface in the far future gives the vacuum persistence 
that explains the Gibbons-Hawking radiation. This raises a fundamental question 
which vacuum persistence amplitude is the proper definition for the effective action 
and the dS radiation. The related issues including the massless scalar field will be 
addressed in a future publication. 
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